Accelerating Expansion of the Universe 
with Nonlinear Spinors 
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Abstract 

The accelerating expansion of the present Universe is an exciting challenge 
to the standard cosmology, which reflects that our current theories on matter 
' have some incompletion. In this paper, we analyze the cosmological model with 

nonlinear spinor fields source in detail. The results may be able to provide a 
natural explanation for the puzzles of the acceleration and negative pressure. 
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1 Introduction 

The observational data show that the expansion of the present Universe is accelerating [1, 
2, 3, 4]. The fundamental quest in theoretical cosmology is to identify the origin of this 
phenomenon. In principle there are two possibilities: (1) the presence of matter has 
negative pressure which violates the strong energy condition, (2) modified Friedmann 
equations [5]. Numerous researches along this direction come forth recently. The follow- 
ing work investigates the properties of a cosmological model with the dominant spinor 
fields as matter source, which is a further research of the preceding papers[6, 7, 8]. 
Some similar works were once provided in [9, 10, 11, 12], where the authors also con- 
sidered the nonlinear spinor field as matter model, but only one spinor field was taken 
into account. 
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In this paper, we translate the dynamic variables of the spnior fields into the equiv- 
alent concepts of mass density p m and pressure P of perfect fluid model, and express 
these variables as well as the other geometrical variables such as the expanding speed 
d, the acceleration d, the Hubble's parameter H, as the functions of scale factor a, and 
then give detailed numerical computation. The influences of the parameters of this 
model on the behavior of the Universe are also considered. The results may be able to 
provide a natural explanation for the puzzles of acceleration and negative pressure. 



2 Dynamic equation 

For the Friedman-Robertson- Walker(FRW) metric 

9ilv = a 2 (t)diag[l, -1, -/ 2 (r), -/ 2 (r) sin 2 0], (2.1) 

with 

{sin r if K = 1 , 
r if K = 0, (2.2) 
sinh r if K = — 1 , 
we derived the dynamic equation with nonlinear spinor source as follows [6] 

a' 2 + Ka 2 - W = ^ E (Kfc + m h ) + C , (2.3) 

where a' = J^a, dVt = f 2 sin 6drd9dip stands for the angular volume independent of t, 
Co is a negative integral constant closely relevant to P, and is the kinetic energy 
of k-th spinor. In statistical sense we have the mean kinetic energy as follows [13] 

K ^ kT( 3mfe 2 

with fh being the mean mass of all particles, and b is a constant determined by the 
initial data 



b 



5kT f 1 + 5kTA <<x (W = 2 .32xl0 4 K). (2.5) 



ao \ m \ Am 

Detail analysis shows we have solution of physical meanings only if K — 1, i.e. the 
closed space. 
Denote 

R = ^FT E m k, C = -R 2 {1 - 5 2 ), (0 < 5 < 1), (2.6) 
611 x k &n 

where R is a cosmic scale independent of a, 5 is a dimensionless number relevant to 
pressure P. Substituting (2.4) and (2.6) into (2.3) we finally get 

a' 2 = \R{sftf + V - a) + 5 2 R 2 + ^Aa 4 - (a - i?) 2 , (2.7) 

o 



we have singular-free cosmological model if 

b<\R{l-5 2 ). (2.8) 
o 

For the extreme case with b = A = 0, we have the simple harmonic solution 

a = R(l - dcost), / = sinr. (2.9) 
Compared with observational data, a rough computation shows 

5 = 0.8-1, R= (0.55 ~ 1.85) x 10 12 yr, t a = 20° ~ 30° (2.10) 
where t a is the present time. 

3 Functions and Properties 

Theoretically, by (2.7) together with the concrete values of constants 5, b, A, we can 
exactly determine the solution. But we will encounter with difficulties in practical 
calculation, because the solution of (2.7) can not expressed by elementary functions 
and the definite values of 5, b, A are unknown. However the problem can be much 
simplified if we take a but not t as independent variable. In what follows, we adopt 
this method and make numerical computation for some important variables. 

We take R — 1 as the length unit. Let A = |Ai? 4 , we get the dimensionless form of 
Eq.(2.7) as follows 

6 



a 



r2 



-{Va 2 + b 2 -a) + 5 2 + Aa 4 - (a - l) 2 , (3.1) 
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5 V« + b 2 

Let r be the cosmic time, i.e. dr = adt, then we have the following dimensionless 
expanding speed a and acceleration d expressed by a 



da 


da 


dr 


adt 


d 2 a 


1 


dr 2 ~ 





± I /6 

a 



(^(Va T Tb 2 - a) + 5 2 + Xa 4 - (a - 1) 2 Y , (3.3) 
a + Aa 4 + 1 - 5 2 ] . (3.4) 



3(^+2^) _2 
5 V^TF 5 

By (3.4) we have 

a — * — ( 1 — S 2 b) — > +oo, as a — > 0. 

a * V 5 / 

So the Universe indeed has a period of inflation, the following numerical results clearly 
show this period. By (2.8), we find the condition of singular-free is just the condition 
of inflation. 
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Now we firstly translate the dynamic variables of the nonlinear spinor model (2.7) 
or (3.1) into the equivalent p m ,P,w of the usual fluid model via the Einstein equation, 
with which all scientists are more familiar, and then look for the relations among them 
and the variable a as well as the parameters 5, b, A. By calculation we have the scalar 
curvature 1Z and Einstein's tensor as 

TZ = 6^^, G 00 = -3(^ + k), Gii = 2--^ + K, (3.5) 
a 6 \a z J a a 1 

and G22 = G11/ 2 , C33 = Gii/ 2 sin 2 #. Substituting (3.5) into Einstein equation with 
fluid model 

T " = ~^ G " + A #) = + - P6 ^ ( 3 - 6 ) 

07TCx 

for static fluid, we get 

Pm = -7^(G° + A) = ^r[3(a' 2 + Ka 2 ) - Aa% (3.7) 
P = 2^h(Gg + n) + 3A] = ^j(a' 2 - 2aa" - tfa 2 + Aa 4 ). (3.8) 

Define the density unit po = s-kGR 2 = ^> anc ^ substituting the dynamic equation of 
spinor model (3.1) and (3.2) into p m , P, then we get the dimensionless functions of the 
nonlinear spinor fields corresponding to the fluid model 



p m = ^ (S 2 - 1 + -(3v / ^T^+ 2a) J , (3.9) 

P=±(p-l + *-£ ), (3.10) 
a 4 V 5^T¥J y 1 



w= p 5V^+¥(5 2 - 1) + Qb 3 

p m 3[6(a 2 + 6 2 ) + % /a^TF(5(5 2 -l)+4a)]' l ' ' 

By the above equation, we find p m , P,w are all continuous functions of a and constants 
5, b, but independent of A as they should be. The Hubble's parameter is defined by 

da a' 1/6 



*=^=?=?u ( ^+^- a)+ * 2+A ° 4 - ( <"- i)2 j • (3 - 12) 

In addition, the angular time t and the cosmic time r can be expressed by the 
integrals of a 

r a 1 r a a 

t = —da, t = / —da, (3.13) 

where a m i n stands for the minimum of a. tH is also an important dimensionless 
variables. According to the functions defined above, we can easily compute the values, 
and clearly describe the properties of the corresponding cosmological model. 

The numerical results are shown in following figures. Fig.(l), (2) and (3) display 
the phase trajectories of a' ~ a with the change of parameters b, A and S respectively. 
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Trajectories of a'~a with varying b 
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Trajectories of a'~a with varying \ 
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Figure 1: with varying b 



Trajectories of a'~a with varying 8 




Dimensionless scale factor a 



Figure 2: with varying A 



Comparison between cosmic time x with angular time t 




Figure 3: with varying 5 Figure 4: Evolution curve & tH 

By these figures we find that, b mainly influences the behavior as a — > a m i n , in contrast, 
A influences the behavior as a — > a max . The influences of 5 on a are more symmetrical. 

By Fig. (4) we find the curve of a(t) is quite near the simple harmonic function. 
The cosmic time r increases slowly as a — > a m i n , but more fast as a — > a max . Except 
for the two extreme cases of a, we find the dimensionless function tH almost acts like 
a line for the most cases. The typical value is 0.5 < tH < 0.7. 

Fig. (5) and Fig. (6) show the Inflating and accelerating expanding process of the 
Universe near a — > a m i n . In contrast, the cosmic time r evolves very slowly in this 
period. These two factors mean that, the Universe just spends a very short time to 
complete her rebirth. 

By Fig. (5) we find H varying slowly except for the case near a m i n . This is the 
reason why the Hubble's law is quite effective. 

Fig. (7) shows the evolution of the equivalent mass density p m and pressure P of 
the spinor fields. In contrast with the fluid model, for the spinor fields we have two 
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Inflation of the Universe 




Dimensionless scale factor a 



Figure 5: The acceleration a(r) & the Hubble's parameter H 

distinct differences: (1) the pressure is negative, (2) the density p m is not a monotone 
decreasing function of a. This means that, a spinor can not be simply treated as a dead 
point particle while the spacetime is heavily curved, which strongly interplays with the 
spacetime for this case. However it becomes normal as the spacetime become flat, e.g. 
as a > 0.2 shown by Fig. (7). 

The curve of the equation of state w — — is shown in Fig. (6), which seems to be 

pm 

quite coincident with the observational data. 

4 Conclusion 

Now we extract some important information from the above analysis and figures: 
(CI). The cosmological factor A have little influence on the young Universe, and 

by Fig. (2) we find that A is purely superfluous without any effective use but increasing 

theoretical troubles, so setting A = is a wise choice. 

(C2). The temperature parameter b increases the singular as a — > 0, but has little 

influence on a normal evolving Universe. Noticing the small value of b, we find (2.9) 

gives a good approximation of the solution. 

(C3). The Universe of nonlinear spinors has an inflating and accelerating expanding 

period in her early age. The inflation means the Universe only spends a very short 

time to get her rebirth and recover her vitality. 
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Figure 7: The relation between p m , P and the scale factor a 
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(C4). Different from the ffuid model, the spinor fields provide negative pressure, 
and the equivalent density p m of the spinor fields is not a monotone decreasing function 
of a. Spinor strongly interplays with the heavily curved spacetime. 
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